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Abstract 

The leading order correction to the metric of a Schwarzschild black hole, due 

to the backreaction of infalling fermionic matter fields, is shown to produce a 

shift of the event horizon such that particles that would constitute Hawking 

radiation at late retarded times are now trapped. Fermionic field operators 

associated with infalling and outgoing modes at the horizon behave canoni- 

eally in the semiclassical approximation. They are, however, shown to satisfy 

a nontrivial exchange algebra given in terms of the backreaction, when the 

shift is 'quantized' by means of correspondence. The consequent exchange 

algebra for bilinear fermionic densities is also obtained. 
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I. INTRODUCTION 



It is well-known that fermionic matter behaves quite differently from scalar quantum 
fields in black hole spacetimes. The foremost difference stems from the fact that fermionic 
fields are not subject to the weak energy condition |T[ T llv ^^ v > on the energy-momentum 
tensor T^ u , where £ A are time-like Killing vectors. An important consequence of this differ- 
ence is the absence of super-radiant scattering of fermions off charged or rotating black holes 
0. The spin statistics theorem which determines the property of the energy momentum 
tensor stated earlier, implies that there can be only one particle per each energy state of the 
fermion. It follows that the amplitude of the reflected wave cannot be enhanced by the black 
hole, thereby preventing the black hole from losing angular momentum or charge by fermion 
emission. The dissimilarity crops up again in the thermal radiation from evaporating black 
holes ||, where the spectral distribution from such a black hole, upon identifying its surface 
gravity with the Hawking temperature, is a Fermi-Dirac distribution for massless fermionic 
fields rather than the Planckian distribution expected for scalars. 

An aspect of Hawking's seminal analysis of black hole radiation that has since received 
attention from diverse angles is the issue of backreaction of matter fields on the black hole 
geometry. In the original (asymptotic) analysis, this was considered only in an adiabatic ap- 
proximation whereby the thermal radiation caused a reduction in the mass of the radiating 
black hole resulting in an increase of the Hawking temperature (and hence in the intensity 
of radiation). The consequent shrinking in the horizon size eventually leads to a complete 
evaporation of the black hole, leaving behind only thermal radiation, - a process which, if 
true, implies a breakdown of standard quantum mechanics with its law of unitary evolution. 
This is the notorious Information Loss Paradox. The resolution of the paradox may in prin- 
ciple entail modification of the principles of quantum mechanics 0, or indeed consideration 
of remnant matter that retains the information (in terms of pure quantum states) that fell 



into the black hole |K|, [11]. A less exotic approach, one which endorses standard quantum 
mechanics and seeks to discern aspects of the original analysis that may be subject to mod- 
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ification upon inclusion of effects neglected earlier, may be another engaging possibility 0] 
- ||. This more conservative viewpoint will be adopted in what follows. 

According to this premise, quantum states characterizing the horizon differ sharply from 
the benign vacuum configuration corresponding to a fixed classical background considered 
in the purely semiclassical approach. The horizon is more likely a dense gas of strongly 
interacting particles at very high energies @], - possibly even a gas of open strings 
||. Infalling particles interact with this system according to the laws of unitary evolution, 
leading to a (possible) description in terms of a quantum mechanical S-matrix Q involving 
only gravitational interactions. The large blue-shifts, close to the horizon, on the solutions 
of the matter field equations in the black hole geometry (which incidentally is germane to 
Hawking's asymptotic analysis H) stand to interpretation in terms of the corresponding 
field quanta acquiring enormous centre-of mass energies. Momentum transfers, on the other 
hand, are taken to be low, on account of the large angular separation expected between the 
most frequently colliding particles. Thus, gravitational interactions tend to dominate over 
all others in these kinematics.^ 

These interactions have been considered for scalar matter fields by Kiem et. al. ||. 
The authors consider linear perturbations of the metric due to matter backreacting classi- 
cally on the geometry. The modified geometry turns out to be the 'cut and paste' of two 
Schwarzschild geometries along the horizon (given by the Kruskal coordinate x + = 0) with 



a shift 6x~ which is given in terms of the classical backreaction ||12|| . This shift, in turn, 
can be transcribed into a shift of the reference null Eddington-Finkelstein (EF) coordinate 
(defining the optimum 'time' at which a massless particle must leave X~ so as to graze the 
horizon and just make it to thereby affecting outgoing Hawking particles: the sign of 
the shift is such that these particles are now trapped. Consequently, the infalling and outgo- 
ing modes of the matter field operator at the horizon, which, in the absence of backreaction 



Admittedly, such kinematical arguments carry a hint of fiat space intuition. 
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are independently measurable observables, are now no longer so. Instead, their commutator 
is proportional to the shift of the reference null EF coordinate. A more nontrivial exchange 
algebra is derived between these modes, by appealing to a correspondence principle to go 
beyond the semiclassical approximation, by promoting the shift to an operator (albeit a 
bilocal one).0 

In this paper, we address the issue of backreaction of neutral massless fermionic matter 
falling on the black hole geometry. In particular, we focus on the anticommutator of the 
infalling and outgoing modes of the fermionic field operator at the horizon of a Schwarzschild 
black hole, to determine whether these modes interact gravitationally just like scalar matter. 
We show that the horizon undergoes a shift quite similar to the scalar case, in a linearized 
Einstein-Cartan formulation. Consistent propagation of fermions in the modified geometry 
leads to a canonical anticommutation relation between these modes in the semiclassical 
approximation. This leads to a canonical algebra for fermionic bilinear densities. Appealing 
next to the correspondence principle, a non-trivial anticommutation exchange algebra is 
derived for the fermionic modes. Fermionic bilinear densities satisfy an exchange algebra 
similar to the one found for scalar matter. 

The paper is organized as follows: In the second section we briefly review the Hawking 
process for fermions. In section 3 we use a perturbed metric to model the effect of the back 
reaction of the fermion energy momentum tensor. In section 4 the effect of this is studied 
and the non-trivial effect on promoting the perturbation as a operator-valued quantum 
fluctuation is investigated. Finally in the concluding section we discuss the limitations of 
the method used, and the issues yet to be addressed in this approach. 



2 Such an algebra between two observables which are usually taken to be mutually commuting 
is strongly reminiscent of Complementarity in quantum mechanics. Hence the name Black Hole 
Complementarity for this result. 
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II. WEYL FERMIONS IN A SCHWARZSCHILD BLACK HOLE 



We begin with the Weyl equations for massless fermions in a Schwarzschild spacetime; 
the metric of the latter is given, in Kruskal-Szekeres coordinates, by 



ds 2 = y 1 e- r ' 2GM dx + dx- - r W . (1) 
r 



?2 = 32(GMf e _ r/2GMdx+dx _ _ ^ w 
r 

A choice of tetrad components e m 01/3 reproducing this metric is given by M 

X+++ x+ (2) 

e ~~ 2GMAV2 y^J 



e x 



2GMAV2 (3) 

e e - + = -1/r e e +- = -l/r (4) 
" + = -z/rsinfl + ~ = z/rsin0 . (5) 

Here, A = (1-2GM). The dotted-undotted pair of indices indicate a tangent space vector 
as usual, with dotted (undotted) indices per se indicating chiral (antichiral) spinors of the 
tangent space Lorentz group. The Weyl equation is then given by, 

iV a/ fya = = *V°% , (6) 

where 

*V a/ Va = ie m > a$ ( d m ^ a - uj a ^) etc. (7) 

and io^ na are the (chiral) spin connection matrices, given in terms of the tetrad components 
by the standard formula 

Near the horizon, A — > ; keeping in mind that length (and time) as measured by asymptotic 
observers scale by the singular factor A in the horizon region, because of infinitely large blue- 
shifts that the solutions of the dynamical equations undergo, one can rescale ip — > ip/A?. 
This reduces the pair (El) to 



d v ^ + = = d v $+ (9) 
d u ^ = = d u ^, (10) 

where u, v are the Eddington-Finkelstein null coordinates. Thus, the Weyl field decomposes 
into 'retarded' (outgoing) and 'advanced' (incoming) solutions near the horizon, similar to 
the scalar field || 

r ut = ^ + («,n) (ii) 
f = ^_(„,n) . (12) 

Since the out (in) solutions constitute a complete set of solutions to the Weyl equation on 
X + (X - ), one can match the advanced and retarded propagation of these solutions at the 
horizon a la' Hawking ||, leading to the reparametrization 

tfj out {u{v),Q) = ip in (v,n) , (13) 

where, u(v) = v — 4GMlog(|^j). The reparametrization is singular at v which represents 
the latest reference 'time' at which an incoming wave leaves X~ to get scattered to X + along 
the event horizon given by u — > oo. For v > Vq, all incoming waves are trapped by the 
black hole. 

Following the original treatment , the reparametrization above is next used to compute 
the (asymptotic) Bogoliubov coefficients appearing in the Bogoliubov transformations con- 
necting the creation-annihilation operators of fields having support on the two asymptotic 
null infinities 

<ko' ( cwtV + Puu'bLui ) , (14) 

where c , c> (a , a* ) are the creation-annihilation operators associated with ip out ( ip in ) on 
X + (X~). The distinction from the scalar case now manifests as one attempts to calculate the 
spectral distribution of the outgoing radiation by calculating the expectation value of the 
number operator c^c in the vacuum on X~; here one must recall that these operators obey 
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an anticommutation algebra instead of a commutation relation. An immediate consequence 
is a thermal spectrum with a Fermi-Dirac distribution, upon identifying the surface gravity 
with the Hawking temperature |§, @. 

This analysis has obviously ignored the change in the black hole geometry induced by 
infalling and outgoing fermionic matter - the backreaction. In the next section we attempt 
to incorporate it in the semiclassical approximation: the linearized change in the black hole 
tetrad components due to infalling fermionic matter is determined, and re-expressed as a 
shift in the horizon. This is then used to determine an exchange algebra for the fermionic 
fields. 



III. CLASSICAL BACKREACTION FOR FERMIONS 

The dominant effect of backreaction |13[of quantum fermionic matter on the classical 
black hole geometry can be characterized by a linearized perturbation of the frame compo- 
nents: e£f — ► e£f + h%P. The linearized fluctuations are related to linearized fluctuations 
of the Schwarzschild metric according to 

where the e£f are the Schwarzschild tetrad components given in (P). We are especially 
interested in the effect of infalling fermionic fields on the black hole geometry. Following 
Hawking's approach to black hole radiation 0, it is sufficient to restrict to waves of very 
high frequency near the horizon, i.e., adopt the geometrical optics approximation. For 
fermion fields falling on the horizon, therefore, the largest contribution to the backreaction 
will come from T x + and T x - the nonzero components of the energy momentum tensor 
in the longitudinal direction, where T%f may be taken to be (the expectation value of) the 
fermionic energy momentum tensor near the horizon. The other energy-momentum density 
components are negligible in the kinematical regime of interest. 
Appealing now to the Einstein-Cartan equation 
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R a J ~ \e a jR = Trf , (16) 

the tetrad fluctuations are seen to obey the linearized Einstein-Cartan equation (see Ap- 
pendix A) 

( V* - 1 )h x+M = k T x+i++ , (17) 
(k is a constant depending upon the mass of the black hole), with the solution 

= J d 2 n' /(fl , fi') T x+i++ , (18) 



where, f(Q,Q') is the Green's function of the Laplacian on the two-sphere | 12| . With this 
modification of the geometry near the black hole horizon, consistent propagation of Weyl 
fermion fields, very close to the horizon, is described by the modified equations (cf. eq.s (|| 
0)) 

aTVs-V- = = x + V x+ ^ + , (19) 
where, V x + = d x + — h x + x +d x - etc. Where h x + x + = et+ h x + , ; . The formal solution for 



ip out = ip + may be written as 

^ + = ^ + ( x - + / dy + h x+x+ (y + ,n)) . (20) 



A similar solution exists for ip m = ip_ with a shift in the other Kruskal coordinate x + . The 
important point to note in (|20|) above is that the integration limit excludes the interval 
(0,Xo ); as pointed out in ||, this region is not interesting for our purpose of estimating the 
effects of backreaction, since for v < Vq, all infalling waves reflect back onto X + . Backreaction 
effects are important only for particles that get trapped behind the horizon. 

Relating the shift 5x~ given in (^) to the affine parameter A of the null geodesic generator 
of the modified event horizon (obtained to linear order in the shift of the tetrad components), 
one can obtain a la' Hawking || a matching condition for v > v o between the incoming and 
outgoing solutions 



S 



iP m (v(u) + 6v ) = V u \u) , (21) 

where, 

v (u) =vo- AGMe (v °- u)/4GM (22) 

5v = - l°° dv [dtt /(fi, n')e {v °- v)/4GM T vv v>v (23) 

Thus, the shift in the optimum value of the EF coordinate on X~ is such that fermions 
that would have made it to X + in the absence of backreaction, are now trapped behind the 
shifted horizon. A remark on the similarity of the shift 5vo found here and that found in 
H for scalar matter is perhaps in order. Recall that near the horizon infalling particles 
are blue-shifted to enormously high energies, so that one can appeal to the geometrical 
optics approximation in dealing with these. Now in flat spacetime, it has been shown |TJj 
that in this approximation (which corresponds to the eikonal approximation), fermionic and 
scalar cross sections become identical for electromagnetic interactions, pointing to an on- 
shell induced supersymmetry. The similarity between the gravitational effect of fermions 
and scalars on the horizon seen here may be attributed to such a supersymmetry. Clearly, 
in the geometrical optics approximation, helicity-flip amplitudes for Weyl fermions vanish. 



IV. EFFECTS OF BACKREATION 

A. Semiclassical Approximation 

The semiclassical approximation, as always, considers quantum (fermionic) matter 
in a classical gravitational background, which in this case is a (backreaction-modified) 
Schwarzschild geometry. In the absence of backreaction, scalar field operators corresponding 
to out and in solutions of the covariant Klein-Gordon equation are known to be mutually 
commuting for v > Vq 0. It is easy to see that, with the modification discussed above to 
the geometry, this is no longer the case: to linear order in the (c-number) shift Svq, one can 
show that 
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[<f) out (u(v), n) , = 2m8v 5(v - v') 5 (2) (^ - fi') . (24) 

Here use has been made of the canonical commutation relation for scalar fields || 

[(f> in (v 1 ,n 1 ) , d V2( p m (v 2 ,Q 2 )} = 27rid{v 12 ) 6® fa -Q 2 ) . (25) 

A similar analysis for fermionic field operators, using the canonical anizcommuting rela- 
tions 

{ i> a (vi, Qi) , ifjpfo, ^2) } = , a,(3 = +,- 

{ , VtK^) } = 27r**(ui 2 ) 5 (2) (^i - fi 2 ) (26) 

leads to a null result to linear order in Svo'. 

{ ^ out , ^ n } = (27) 

as one might have expected. Indeed, fermionic fields are themselves unobservable; densities 
constructed out of fermionic bilinears are the true observable quantities. The simplest 
bilinear composites constructed out of fermions are the components of the current in the 
tangent plane J a a. The dominant components of the current are0 

J m ( V ) EE EE 1>_{V)^{V) , J° Ut (u) EE ^+ , (28) 

These bilinear observables turn out also to behave canonically, i.e., [J out , J m ] = to lowest 
order in 5vq (see Appendix B for a sketch of the commutator calculation). This behaviour 
is quite in contrast to the behaviour of scalar fields in the semiclassical approximation: the 
infalling and outgoing modes of the fermion field appear to be completely independent, i.e., 
non-interacting, despite the classical shift of the horizon. 

3 The other components are negligible in the kinematical situation under consideration. Also, 
scalar bilinears vanish for chiral fermions. 



10 



B. Quantum Correspondence and Exchange Algebras 

With fields as operators, the energy-momentum flux operator^ P v = dv T vv has a 
nontrivial commutation relation with the incoming field. It generates translations along the 
v direction. 

[p v {v), i? n {v')\ = 2mS 2 (n - n')d v i} in (v', n') (29) 

Now, as seen in equation (23), the shift 5vq is related to the energy momentum tensor. If we 



assume that this classical relation can be promoted as such to a relation between operators 
by means of some sort of a quantum correspondence principle 0, it can be used find the 
commutation relation of Svq with the incoming field 

[6vo{n) , ip in (tt')] = -167iif(n,n')e {V0 - v)/4GM d v ^ in v > v . (30) 

This result can now be used, following Kiem et. al. |J to determine an exchange algebra 
between the in and out field operators. Keeping in mind the canonical anticommutation 
relations obeyed by the fermion fields, this algebra (to lowest order in the backreaction) is 
given by 

ip out (u,n)?p m (v,n') = -{ i - i6mf(n,n')e {u - v)/4GM d u d v })i} in (v,n')?p out (u,n) . (31) 

In the absence of backreaction f(Q, Q') = 0, we get the standard anticommutation relation 
between the field operators. 

Thus to the extent one can trust the procedure of promoting ([23|) to the level of an 
operator relation, fermionic field operators appear to obey the requirements of Complemen- 
tarity. The point becomes clearer when one computes, within the same correspondence 
approach, the exchange algebra of fermionic bilinear densities; this turns out to be similar 
to the exchange algebra of scalar fields 



J The Energy Momentum tensor is not well-defined near the horizon, but the momentum flux is 



11 



j out (u,n)j in (v,n') = { 1 - iQmf(n,n')e iu - v)/iGM d u d v }j in (v,n')j out ( u ,n) . (32) 

In other words, observables, mutually commuting in absence of backreaction, now obey a 
nontrivial exchange algebra (for v > vq). The similarity of the exchange algebras for the 
fermion bilinears and scalar fields is of course, no surprise, being expected on account of the 
similarity in the behaviour of fermions and scalars under spacetime translations which lies 
at the root of the derivation of these algebras. In addition, as already mentioned, since the 
geometrical optics approximation is essentially being used here, and fermions and scalars 
behave similarly in this approximation, the only difference that may be expected is from 
the statistics, i.e., the fact that fermionic field operators obey a canonical anticommutator 
algebra. This is as obvious here as in the case of the spectral distribution of Hawking 
radiation. 

V. CONCLUSIONS 

In the semiclassical approximation, backreaction of matter does not seem to induce any 
interaction between infalling and outgoing modes of the fermion field, while for scalars such 
interactions do indeed appear. It is important to keep in mind however that in this analysis, 
tranverse gravitational interactions as well as non-gravitational forces have been ignored on 
the plea that they would be subdominant in this situation. Such forces could conceivably 
change this results somewhat J7|, although a complete picture can only be expected after 
nonperturbative 'quantum gravity' is understood. In any case, our analysis does not indicate 
any serious flaw with the semiclassical approximation in the fermionic case, even though for 
scalars going beyond this approximation seems imperative. A deeper understanding of this 
difference, possibly in terms of foliations of the geometry (e.g., along lines suggested in |15*| ) 
is certainly desirable. It might be related to the diffrence mentioned at the outset, namely 
that fermions do not obey the weak energy condition. 

As mentioned earlier, the validity of the exchange algebras derived above hinges on the 
assumption that the shift Svq, while a shift in the optimum value of a coordinate, can be 

12 



elevated to the level of an operator on Fock space. Clearly, as an operator, 5vq is bilocal 
which might possibly underlie a justifiable suspicion of a violation of microcausality. In the 
scalar field case, Kiem et. al. introduce a third scalar field <p hor to ensure that microcausality 
is maintained, with <fr tn evolving unitarily to it at the horizon from T~ . It is not clear whether 
this is a valid procedure without a full quantum theory of gravity to back it up. We therefore 
adopt a deliberately ambivalent stand on this issue. 

In fact, even promoting equation (|23"D to the level of an operator equation and thereby 
endowing operatorial status to a coordinate, purely by appealing to a correspondence prin- 
ciple, is questionable in a situation where a well-defined quantum theory is yet unknown. 
Usually, a correspondence principle works best in the reverse manner, i.e., in discerning the 
structure of a classical theory corresponding to a known quantum theory. The application 
of correspondence in the present context (which elevates the linearized Einstein equation 
to the level of an operator equation), while going beyond the purely semiclassical approx- 
imation (where gravity is classical), appears to be replete with ambiguities, again in the 
absence of a proper 'quantum geometry'. On the other hand, the super-Planckian energies 
acquired by field quanta near the horizon does seem to warrant an approach beyond the 
semiclassical approximation. In this context, a better approach than the one adopted here 
might be to follow the quantization of monopoles in gauge field theory, using the method 
of collective coordinates. Indeed, the horizon, as already mentioned, is modeled by a dense 
two dimensional gas of strongly interacting particles. Such a system, like a large nucleus, is 
very likely to undergo collective motion which should be amenable to exact analysis, unlike 
the (linearized) metric (or tetrad) fluctuations dealt with perturbatively here (the latter are 
analogous to Gaussian fluctuations around the ground state in the case of the nucleus) . We 
hope to report on these issues elsewhere. 
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APPENDIX A: 



The linear order perturbation to the tetrads is shown to satisfy equation (17) here. The 
tetrad is perturbed to linear order e£f — ► e^f + . Hence, to linear order in h, the two 
form, 

R a tt = dm Surf'* - d n 8utf<* + 8utf> + d$> a '^S^ af$ , 

- Su^ uf' a '^ - Suf > a ' $ ' . (Al) 

m,a'p' n m,a'fJ' 71 v ' 

Where 5u stands for the linearised spin connection. The change in the tetrad, is only along 
the longitudinal direction as near the horizon the energy momentum tensor of the matter 
field is dominant in these directions only. 

Se x+t+ + = h x+)+ +(x + ,x-,n), 5e x -_^ = h x ,_^(x + ,x-,n) (A2) 

The changes in the spin connection very near the horizon (neglecting terms proportional to 
x + and higher powers), are calculated and the relevant self dual parts are: 

(^)± ± = (A3) 
(5u x+ )- = -- (do + h x+t++ (A4) 

Similar equations are obtained for 5u x - . The spin connections of the original metric were, 
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(o; a .+) ± 



1 4(GM)' 
2 X + r 2 



1 4(GM) : 
2x~ r 2 



(A5) 



± T = =F A 1/2 , w^jf = T «cos0, ^ >± T = i sin 6>A 1/2 . 



(A6) 



Using the above the perturbed curvature can be calculated from equation Al. The contrac 



tion with the vierbein gives -R m7< 5 
equation fllTD follows. 



van, af3j6 



. From Einstein-Cartan equation ( |iq 



APPENDIX B: 



The commutator of the currents are calculated to order Svq. 



jout jin 



= Vm- i 1 -} V>- - f^+, if>-} ip+tp- 

+ - ^•{^+,'0-}V+- (Bl) 

Using the anticommutators given in equations ( f26|) and the relation of the incoming and 
outgoing waves, 



^ = ^ + ( U ) = ^_{y' + SV ) , $ in = l/i_(v) , 

the above commutator is simplified to: 

2m5 2 (n'-n) [jj+i)- - tfi^+] [5(v ' - v) + 5v d v ,5(v' - v)] . 

Which can be further simplified by using equation( p2|) , to 

~J out , J in ] = 2m5{9! - n)Svo [d v >S(v' - v) (^(v')^-(u) - ^-{v) 
+ 5(v'-v) (d v 4-(v')ip-(v) - ^(v)^(v')J 
= 2m8{Q! - n)8v Q d v , [5(v' - v) (f^(v')^(v) - ip^{v)^{v' 
= . 

Hence in the semiclassical approximation the commutator is the canonical one. 



(B2) 



(B3) 



(B4) 
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